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It is well-known from numerical and experimental results that the fraction of Rydberg states 
(excited neutral atoms) created by tunnel ionization declines dramatically with increasing ellipticity 
of laser light. We present a method to analyze this dependence on ellipticity, deriving a probability 
distribution of Rydberg states that agrees closely with recent experimental [l] and numerical results. 
In particular, contradicting the existing proposed mechanism 1 and general expectations, our results 
indicate that the dependence of Rydberg yield on laser ellipticity is not caused by a rescattering 
process. Rather, we show that most Rydberg electrons never come back to the vicinity of the exit 
point after ionization, and end up relatively far from the atom (compared to the exit point) after the 
laser pulse has passed. We also present experimental data revealing a bifurcation that corresponds 
to a cut-off in Rydberg generation, and present an analytic derivation with a perturbative inclusion 
of the Coloumb force. 
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Tunnel ionization which can occur when an ion- 
izing laser field strength is comparable to the Coulomb 
force binding the outer electron to its atom, is behind 
many recent break-throughs in attosecond science. The 
3-step model, which involves tunnel ionization, classical 
propagation and recollision [3 , works surprisingly well 
in clarifying many phenomena, most notably High Har- 
monic Generation (HHG) [HIS] and Non- Sequential Mul- 
tiple Ionization (NSMI) [6]-[8] . Recently, the role of 
Coulomb force in laser ionization has drawn considerable 
attention from the strong field atomic physics commu- 
nity. In particular, it is believed that to properly study 
many phenomena (e.g. focusing [9 , low energy struc- 
tures at mid-infrared wavelengths [10 , certain asymme- 
tries in momentum spectra [11 , tunneling geometry [TT 
and most notably for our analysis, the generation of Ry- 
dberg states [U |13]), the inclusion of Coulomb force is 
essential. 

The study of Rydberg state creation has attracted re- 
cent interest, with the first direct experimental obser- 
vation of these excited neutrals presented in [1]. The 
main conclusion was said to provide "strong experimental 
support" for a rescattering process, grouping it together 
with HHG and NSMI [T. In contrast, our results using 
only the standard assumptions from strong field physics 
(namely the semiclassical tunneling model [HE]) show 
close quantitative agreement with the same experimental 
data without the need to account for rescattering. More- 
over, in contradiction to expectation that the inclusion 
of Coulomb force on the dynamics is essential, we analyt- 
ically show that Rydberg yield dependence on ellipticity 
can be closely reproduced while neglecting the influence 
of the Coulomb force on the trajectory during the entire 
duration of the long laser pulse. 
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Rydberg states can be created when the tunneled elec- 
tron does not gain sufficient energy from the laser pulse 
and is subsequently captured by the Coulomb field, cre- 
ating an excited neutral atom. While some progress has 
been made in analyzing the generation of Rydberg states 
in linearly polarized light [TSl [15] , the role of polariza- 
tion, and in particular the dramatic decline of Rydberg 
states with increasing ellipticity of light remained to be 
quantifiably explained. Here, we correct and clarify this 
dependence of Rydberg yield on ellipticity, deriving a 
probability distribution which shows excellent agreement 
with the above mentioned experimental [1 and numerical 
data, and which furthermore quantifies the dependence 
of Rydberg states on laser field intensity. 

We then present our own experimental data showing 
a hitherto unobserved bifurcation structure in electron 
energy distribution at the detector. This bifurcation in 
the electron energy from zero coincides with the cut-off 
for generation of Rydberg states for electrons born in the 
continuum with zero velocity (corresponding to the peak 
of the probability distribution). We show that the value 
of ellipticity at which this bifurcation occurs can be an- 
alytically derived and depends on the interplay between 
the Coulomb force and the laser field. 

The electric field of a laser propagating in the z- 
direction is given by: 

F(t) = [cos {ut + 0c) ^ + esin {ujt + 0c) y] (1) 

Ve^ + 1 

where uj is the frequency of the laser, e is the ellipticity 
(the major axis of polarization is along f), and f{t) is 
the slowly- varying pulse envelope: fmax = /(O) = 1- 
For many-cycle pulses, we can approximate the carrier- 
envelope offset (CEO) phase [16]: 0c ^ 0- 

To lowest order, the electron dynamics following tun- 
nel ionization are given in atomic units by [3]: d^f/dt^ = 
—F{t). This simple equation is based on the assumption 
within the Strong Field Approximation (SFA) [S] [17], 
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that the Coulomb force after ionization is neghgible com- 
pared to the laser field: |F| ^ QIt. For the evalua- 
tion of Rydberg states, which by definition depend on 
the Coulomb potential, we will therefore only neglect the 
Coulomb force up till the end of the laser pulse. After 
the laser pulse has passed, for t > Tp (where Tp is the 
ulse duration), we have f{t) ~ 0, so that from Eqn. 
jj): F{t) ~ and therefore the Coulomb force is no 
longer negligible compared to the laser field, capturing 
lower energy trajectories into Rydberg states. 

For a slowly- varying envelope, where df(t)/dt <C cj, 
the force equation can be integrated to give approximate 
dynamics for electrons following ionization, 

/ N ^0 (1 - f(t)cos(ujt)) 



yield (arb. units) 



yit) 
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where At = t — ti (with ti being the instance of ioniza- 
tion) , and Xe is the exit point from the tunnel P [12] . As 
an initial condition, we used Vx (U) = (in accordance 
with the tunneling model [14 ), which is valid for small 
e. The parallel velocity after the laser pulse has passed, 
Vxo, is determined by the phase of the field, at the 
instance of ionization: Vxo = — ^iFo/cjVl + e^. 

The transverse velocity at the exit point is given by 
('^7/5 ^z) and has a Gaussian probability distribution given 
by the well-known ADK formula [14], with standard de- 
viation (SD): 



(4) 



where 7 is the Keldysh parameter: 7 = ujy^2Ip/F. 

Equation Q shows that an electron that tunnels near 
the peak, where Vxo ^ 0, will continue, when averaged 
over a single cycle, to get farther from the exit point dur- 
ing the duration of the pulse. This is because f{t) is a 
monotonically decreasing function for t > 0. In partic- 
ular, this means that most Rydberg states, which form 
predominantly near the peak of the laser field (with for- 
mation at the peak corresponding to a dominant quan- 
tum number of n=8 for Helium [1]), will not rescatter. 

The overall structure of the following derivation for 
Rydberg yield is: We show that the electrons end up 
far from the exit point, Xg, after the laser pulse has 
passed. This imposes a near zero final momentum condi- 
tion. While the momentum along x is close to zero due to 
ionization near the peak, the final transverse momentum 
at the exit point needs to cancel out the drift velocity 
created by the laser field. Therefore, the probability of a 
Rydberg state can be obtained from the probability that 
the transverse velocity at the exit point (which is given 
by the ADK formula [14 ) is such that it approximately 
cancels the deterministic drift velocity. 

From Eqn. ([2|, the position of the electron ionized 




FIG. 1: Analytic curve (red line), simulations (•) and 
experimental data (A) for the Rydberg yield. Total ion 
yield of i^e+, given by Pq, is shown on the same plot 

(blue line) . Experimental data is taken from [1 . 

Tp = 30/s, uj = 0.056 a.u.. Intensity = IPW/cm'^ 



near the peak after the pulse has passed is: 



Xf 



(5) 



The right-hand side of the equation holds for small e, and 
in the framework of the strong field infrared laser physics 
where Fq/u'^ Xq. 

From Eqn. (|5|, it is clear that, in contrast to ionization 
by a plane wave, an electron ionized near the peak of a 
laser pulse will end up a significant distance from the 
parent atom after the laser pulse has passed. Since the 
distance from the parent atom at the end of the pulse is 
large, Xf ^ Xg, the velocity has to be small to satisfy 
the negative energy condition: —Q/xf + KE < E < 0, 
necessary for Rydberg states. We therefore require that 
the transverse velocity of the electron after the pulse has 
passed satisfies the following condition. 



Vy =Vy^ Vdrift ^ 

where Vy was obtained from Eqn. ([3|, and 

Vdrift 



eFo/uj\/^- 



1 



(6) 



(7) 



is the drift velocity caused by the vector potential of the 
laser field [3 [m 

From Eqn. Q, the scaled probability of generating 
an electron with Vy ^ —Vdrifti satisfying Eqn. (|6|, 
is PR{e) ^ Poexp {-e'^Fo{2Ip)^^'^ /u'^) , where Pq is the 
ionization probability in the tunneling regime along the 
major axis of polarization [2 . In the range of small 
it can be approximated to exponential accuracy as: 
= exp (-2(2/p)3/Vl + e73Po) [18- The exponent 



Po 

of Po can be further Taylor expanded to obtain. 



Puie) 
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where is the standard deviation of the Gaussian prob- 
ability distribution for Rydberg state yield as a function 
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FIG. 2: a) Energy of electrons observed at the detector 
as a function of e. b) Corresponding momentum, py 
along the minor axis of polarization. Tp = 33/s, 
Ao 788nm, Intensity {).^PW/cm^. 




FIG. 3: (a) Monte Carlo simulations showing a 
bifurcation of Py near |e| ^0.1, in agreement with 
experimental results in Fig. [2](b). (b) Bifurcation plot. 
Red: Coulomb correction to the momentum, pcyi given 
by Eqn. (10). Blue: drift momentum given by Eqn. 



Bifurcation is estimated to occur at the intersection. 



of ellipticity of light: 



3 + 7' 



(9) 



In deriving Eqn. ([8| a Taylor expansion was used, retain- 
ing only terms up to the order of e^, since <C 1 (the 
Rydberg state generation is only found at small values of 
e). In the tunneling regime where 7 <C 1, Eqn. ([9| can 
be approximated as: <Je = uj / ^/2FQ{2IpY^^ . 

The SD given in Eqn. ([9| increases with decreasing Fq, 
predicting a slower decline of Rydberg states with e, as 
laser intensity goes down. This is in contrast to the total 
ionization yield, given by Pq, where smaller laser inten- 
sities lead to a faster decline with ellipticity. These two 
factors lead to an increase of Rydberg trajectories relative 
to the total ionization yield, explaining and quantifying 
the findings in [1 , where numerics show "an increasing 
percentage of bound trajectories ... with decreasing laser 
intensity." On the other hand, for linearly polarized light, 
where e = 0, Eqn. ([8| predicts an increase in Rydberg 
states with increasing intensity, also in agreement with 
the experimental and numerical findings in [T]. 

Monte Carlo simulations were performed whereby an 
ensemble of 10^ ionization events with a randomly cho- 
sen CEO phase were generated for each value of e over 
the interval [-1,1]. The ionization rate for each event is 



given by Tong et al. [19], assuming instantaneous tun- 
neling delay time [20]. The equation of motion is solved 
with the initial condition that the electron starts its tra- 
jectory outside the potential barrier, at the exit point, 
Xg, with an initial longitudinal momentum of zero and a 
transverse momentum distribution given by Eqn. Q. 

A comparison between the experimental data (taken 
from [1 ), simulations and analytically obtained curve 
given by Eqn. ([8| is shown in Fig. [l] The intensity 
and laser frequency used in the figure were chosen for 
the purposes of comparison with the experiment in 
While the total ion yield varies slowly with increasing el- 
lipticity, the yield of Rydberg states drastically decreases 
to essentially zero for e > 0.3. This holds true over a wide 
span of intensities (0.35 — 3PVK/cm^), where analysis (see 
Eqn. (|8|) and numerics show the disappearance of Ryd- 
berg states in the 0.2 < e < 0.3 range, depending on field 
strength. This dramatic decrease with ellipticity is also 
typical of a rescattering process as measured for HHG 
and NSMI [H |6] . 

Experiment Our experiment provides indirect infor- 
mation on Rydberg state generation by measuring the 
yield of £^ ~ states at the detector as a function of 
e. £^ ^ states can be considered as a cut-off for Ryd- 
berg states since the condition defining the Rydberg state 
(with the corresponding quantum number, n) is given by: 
E = -l/2n^ < 0. 

The experimental setup was as follows: a laser pulse of 
duration Tp = 33/5, central wavelength Aq = 788nm, and 
peak intensity of 0.8PW/cm'^ (CEO phase [16 was not 
stabilized) was produced by a Ti: Sapphire laser system 
focused onto helium atoms in a cold gas jet, with the gas 
jet density adjusted such that on average much less then 
one ionization occurs per laser shot. A COLTRIMS setup 
[20 measures the ion momentum, which is the negative of 
the electron momentum due to momentum conservation. 
The momentum resolution is 0.1 a.u. in time-of-fiight 
direction and 0.9 a.u. in gas jet direction, mainly de- 
termined by thermal spread. A broadband quarter-wave 
plate is used to alter the ellipticity of the laser pulses. In 
the final analysis, the ellipticity and the angular orien- 
tation of the polarization ellipse are calculated from the 
angle of the quarter- wave plate, respecting its wavelength 
dependence [21 . The knowledge of the ellipticity and 
the angle of the polarization ellipse for each detected ion 
allow generating ellipticity-resolved spectra with a high 
resolution. 

The experimental results are presented in Fig. [2j Fig- 
ure [2] (a) shows that E ^ states decline dramatically 
and virtually disappear for e > 0.3, in agreement with 
theory and numerics, shown in Fig. [l] which also predict 
the disappearance of Rydberg states close to e ~ 0.3. 

The experimental data in Fig. [2] (b), and the corre- 
sponding Monte Carlo simulation in Fig. [3] (a), show 
an interesting bifurcation behavior where the center of 
the Py distribution, where Py is the momenta along the 
minor axis of polarization, splits from ^ as e is in- 
creased. (By bifurcation, we mean a qualitative change 
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in behavior as one of the parameters, in this case e, is 
increased.) This bifurcation can not be explained within 
the SFA approximation, which neglects Coulomb force 
after ionization, predicting a close to linear increase in 
\py \ as |e| increases [22 , which would result in immediate 
splitting as e is increased from zero. To explain the bifur- 
cation behavior shown in Fig. |2](b), we therefore have to 
take account of the Coulomb correction to the electron 
dynamics after ionization. Using perturbative approach 
pT] . the total Coulomb correction to the momentum is 
estimated by assuming that to lowest order the electron 
follows a trajectory given by Eqns. ([2| and ([3|. 

The most probable trajectory corresponds to ioniza- 
tion at the peak with v = 0^ and is given by ro{t) = 
xo{t)x -\-yo{t)y^ where xo{t) and ^o(^) are given by Eqns. 
([2| and ([3| with Vx = Vy = 0. The total Coulomb cor- 
rection to the center of the distribution is [11 , 



dti-j^m±M}y]=p^^,+p^^y (10) 



where the integration is from the time of ionization at 
the peak, where = 0, till the end of the pulse at t ^ r^. 

The total momenta along as observed at the detec- 
tor, is: Py = Vdrift-^Pcy A non-zero drift of the center of 
the distribution occurs for \py \ > 0, with the correspond- 
ing bifurcation around: Vdrift = —Pcy Using Eqns. (|3|, 
^ and ([10|, 



- / dt- 
^ Jo ( 



ujt — f{t)sin{ujt) 
{xl{t,eb)^y'^{t,eb)) 



3/2 



(11) 



where, as before, xo(t, 65) and yo{t^ 65) are given by Eqns. 
([2| and ^ with Vx, Vy, and ti set to zero. 

Equation ( pTj ) can be solved numerically for the bifur- 
cation value, 65. Using the same parameters as for the 
experimental data shown in Fig. [2j we get a bifurca- 
tion value of 65 ~ 0.06. This estimate is in approximate 
agreement with experimental data, (Fig. [2](b)) and with 



numerics (Fig. [s] (a)), which both show a bifurcation 
near e ^ 0.1. As shown in Fig. 3 (b), the bifurcation 
point corresponds to the value of e where the Coulomb 
force becomes too weak to offset the velocity drift along 
the minor axis of polarization created by the laser field. 
An under-estimate of 65 is expected when using a per- 
turbative approach since an electron that follows the un- 
perturbed path, fo, feels no Coulomb force and therefore 
gets away from the atom faster than the actual trajec- 
tory, thereby incurring less of a Coulomb correction. 

In conclusion, we analyzed the dependence of Rydberg 
states on laser ellipticity, deriving a Gaussian probability 
distribution for the yield of neutral Rydberg atoms as a 
function of e; this puts recent experimental and numerical 
results into a theoretical framework. In particular, our 
work suggests that rescattering does not play a significant 
role in the creation of excited neutrals, in contradiction 
to the mechanism proposed in [1 . 

Our analysis predicts the disappearance of Rydberg 
states for polarizations with e > 0.3 (when the electron 
yield from tunnel ionization is still quite substantial). 
The standard deviation of the derived probability dis- 
tribution, given by Eqn. ([9]), is found to decrease with 
increasing laser intensity, in agreement with the observed 
decline of Rydberg states at higher laser intensities. 

The experimental data presented here confirms the dis- 
appearance of Rydberg states for e > 0.3 by measuring 
zero energy states at the detector. These states can be 
considered as a cut-off for the creation of bound trajecto- 
ries since they correspond to asymptotically large values 
of n (the quantum number of the Rydberg state). In ad- 
dition, our experiment shows a bifurcation from zero in 
the final momentum, py^ as ellipticity is increased. This 
bifurcation was explained by a perturbative inclusion of 
a Coulomb correction. 
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